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Abstract Although general relativity (GR) passes all present experimental tests with 
flying colors, it remains important to study alternative theories of gravity for several 
theoretical and phenomenological reasons that we recall in these lecture notes. The 
various possible ways of modifying GR are presented, and we notably show that 
the motion of massive bodies may be changed even if one assumes that matter is 
minimally coupled to the metric as in GR. This is illustrated with the particular case 
of scalar-tensor theories of gravity, whose Fokker action is discussed, and we also 
mention the consequences of the no-hair theorem on the motion of black holes. The 
finite size of the bodies modifies their motion with respect to pointlike particles, and 
we give a simple argument showing that the corresponding effects are generically 
much larger in alternative theories than in GR. We also discuss possible modifica- 
tions of Newtonian dynamics (MOND) at large distances, which have been proposed 
to avoid the dark matter hypothesis. We underline that all the previous classes of al- 
ternatives to GR may a priori be used to predict such a phenomenology, but that 
they generically involve several theoretical and experimental difficulties. 



1 Introduction 

Since general relativity (GR) is superbly consistent with all precision experimental 
tests — as we will see below several examples, one may naturally ask the ques- 
tion: Why should we consider alternative theories of gravity? The reason is actually 
threefold. First, it is quite instructive to contrast GR's predictions with those of alter- 
native models in order to understand better which features of the theory have been 
experimentally tested, and what new observations may allow us to test the remaining 
features. Second, theoretical attempts at quantizing gravity or unifying it with other 
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interactions generically predict the existence of partners to the graviton, i.e., extra 
fields contributing to the gravitational force. This is notably the case in all extra- 
dimensional (Kaluza-Klein) theories, where the components g^b of the metric tensor 
[where a and b belong to the D — 4 extra dimensions] behave as (D — 4) (D — 3)/2 
scalar fields (called moduli) in four dimensions. Independently of such moduli, su- 
persymmetry (notably needed in string theory) also implies the existence of several 
fields in the graviton supermultiplet, in particular another scalar called the dilaton. 
The third reason why it remains important to study alternative theories of gravity 
is the existence of several puzzling experimental issues. Cosmological data are no- 
tably consistent with a Universe filled with about 72% of "dark energy" (a fluid with 
negative pressure opposite to its energy density) and 24% of "dark matter" (a fluid 
with negligible pressure and vanishingly small interaction with ordinary matter and 
itself) mOOl llOfl . Another strange phenomenon is the anomalous extra accelera- 
tion towards the Sun that the two Pioneer spacecrafts have undergone beyond 30 
astronomical units P1 l84l[T06l (see also lf73l ). Although such issues do not threaten 
directly GR itself, since they may be explained by the existence of unknown "dark" 
fluids (or by yet unmodelled sources of noise in the case of the Pioneer anomaly), 
they may nevertheless be a hint that something needs to be changed in the gravita- 
tional law at large distances. 

Since many theoretical and experimental physicists devised their own gravity 
models, the field of alternative theories is much too wide for the present lecture 
notes. Detailed reviews may be found in Refs. 11 121 11131 11071 ISTl . We shall focus 
here on the particular case of scalar-tensor theories and some of their generaliza- 
tions. Before introducing them, let us recall that GR is based on two independent 
hypotheses, which can be most conveniently described by writing its action 



where g denotes the determinant of the metric g^v, R its scalar curvature, and we 
use the sign conventions of Ref. fSTI, notably the mostly-plus signature. The first 
assumption of GR is that matter fields are universally and minimally coupled to one 
single metric tensor g^y. This ensures the "Einstein equivalence principle", whose 
consequences will be summarized in Sec. |2] below. The second hypothesis of GR 
is that this metric g/jy propagates as a pure spin-2 field, i.e., that its kinetic term is 
given by the Einstein-Hilbert action. The core of the present lecture notes, Secs.|3]to 
|6] will be devoted to the observational consequences of other possible kinetic terms. 



2 Modifying the matter action 

In the above action (O, square brackets in ^matter [matter, g^jy] mean a functional 
dependence on the fields, i.e., it also depends on their first derivatives. For instance. 
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depends not only on its spacetime position x but also on its 4-velocity = dx^ /dt. 
Since the matter action defines the motion of matter in a given metric g^^, it is a 
priori what needs to be modified with respect to GR in order to predict different 
trajectories. This idea has been studied in depth by Milgrom in 1791 [80l . where he 
assumed that the action of a point particle could also depend on its acceleration 
and even higher time-derivatives: 5'pp(x, v,a,a, . . .). However, any modification with 
respect to action (|2]i is tightly constrained experimentally for usual accelerations, 
notably by high-precision tests of special relativity. On the other hand, physics may 
happen to differ for tiny accelerations, much smaller than the Earth's gravitational 
attraction. In such a case, the mathematical consistency of the theory may be in- 
voked to restrict the space of allowed theories. A theorem derived by Ostrogradski 
in 1 850 ISTlll 161 shows notably that the Hamiltonian is generically unbounded from 
below if 5pp(x, v,a, . . . ,d"x/dt") depends on a finite number of time-derivatives, and 
therefore that the theory is unstable. A possible solution would thus be to consider 
nonlocal theories, depending on a infinite number of time derivatives. This is actu- 
ally what Milgrom found to be necessary in order to recover the Newtonian limit 
and satisfy Galileo invariance. Although nonlocal theories are worth studying, and 
are actually obtained as effective models of string theory, their phenomenology is 
quite difficult to analyze, and we will not consider them any longer in the present 
lecture notes. General discussions and specific models may be found for instance in 
Il98l|79l|99l|49l|4l. 

Another possible modification of the matter action ^matter [matter, g^v] is actually 
predicted by string theory: Different matter fields are coupled to different metric ten- 
sors, and the action takes thus the form S'matter [matter''', ^|jy]. In other words, two 
different bodies a priori do not feel the same geometry, and their accelerations may 
differ both in norm and direction. However, the universality of free fall is extremely 
well tested experimentally, as well as the three other observational consequences of 
a metric coupling S'matter [matter, g^y], that we will recall below. The conclusion is 

that string theory must actually show that the different metrics gj^y almost equal 
to each other One possible reason is that their differences may be mediated by mas- 
sive fields, and would become thus exponentially small at large enough distances. 
But even in presence of massless fields contributing to the difference between the 
various g|jy, a generic mechanism has been shown to attract the theory towards GR 
during the cosmological expansion of the Universe 1441 l46l l47l . 

Let us now recall the four observational consequences of a metric coupling 
5matter[matter,g^v], as well as their best experimental verifications. If all matter 
fields feel the same metric g|,y, it is possible to define a "Fermi coordinate system" 
along any worldline, such that the metric takes the diagonal form diag(— 1, 1,1,1) 
and its first derivatives vanish. In other words, up to small tidal effects proportional 
to the spatial distance to the worldline, everything behaves as in special relativity. 
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This is the mathematically well-defined notion of a freely-falUng elevator. The ef- 
facement of gravity in this coordinate system implies that (i) all coupling constants 
and mass scales of the Standard Model of particle physics are indeed space and time 
independent. One of the best experimental confirmations is the time-independence 
of the fine-structure constant, \oi/a\ < 7 x 10^'^yr^', six orders of magnitude 
smaller than the inverse age of the Universe 1971 [37l l58l . A second consequence 
of the validity of special relativity within the freely-falling elevator is that (ii) lo- 
cal (non-gravitational) experiments are Lorentz invariant. The isotropy of space has 
notably been tested at the 10^^^ level in 1881 l74l l33l . The third consequence of a 
metric coupling is the very existence of this freely falling elevator where gravity 
is effaced, i.e., (iii) the universality of free fall: (non self-gravitating) bodies fall 
with the same acceleration in an external gravitational field. This has been tested 
at a few parts in lO'^^ both in laboratory experiments l9||2l, and by studying the 
relative acceleration of the Earth and the Moon towards the Sun II 151 . The fourth 
consequence of a metric coupling is (iv) the universality of gravitational redshift. It 
may be understood intuitively by invoking the equivalence between the physics in 
a gravitational field and within an accelerated rocket: The classical Doppler effect 
suffices to shows that clocks at the two ends of the rocket do not tick at the same 
rate, and one can immediately deduce that lower clocks are slower in a gravitational 
field. More precisely, a metric coupling implies that in a static Newtonian poten- 
tial goo = — 1 + 2t/(x)/c^ + /c*), the proper times measured by two clocks is 
such that Ti /t2 = 1 + [C/(xi ) - t/(x2)] /c^ + /c'^). This has been tested at the 
2 X 10^^ level thirty years ago by flying a hydrogen maser clock II108II1091 . and the 
planned Pharao/Aces mission ||9T1 should increase the precision by two orders of 
magnitude. 

In conclusion, the four consequences of a metric coupling have been very well 
tested experimentally, notably the universality of free fall (i.e., the relative motion 
of massive bodies in a gravitational field). Therefore, although theoretical consid- 
erations let us expect that the Einstein equivalence principle is violated at a funda- 
mental level, we do know that deviations from GR are beyond present experimental 
accuracy. In the following, we will thus restrict our discussion to theories which 
satisfy exactly this principle, i.e., which assume the matter action takes the form 
^matter [matter, g^v]- On the other hand, we will now assume that the kinetic term of 
the gravitational field, say ^gravity, is not necessarily given by the Einstein-Hilbert 
action of Eq. ([T|i- 



3 Modified motion in metric theories? 

For a given background metric g^^^, the kinetic term ^gravity defines how gravita- 
tional waves propagate, and the matter action S'matter how massive bodies move in 
spacetime. If we assume a universal metric coupling ^matter [matter, g^v] as in GR, 
we are thus tempted to conclude that the motion of matter must be strictly the same 
as in GR, and that the present lecture notes should stop here. However, S'gravity also 
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defines how is generated by the matter distribution. Therefore, the motion of 
massive bodies within this metric does actually depend directly on the dynamics of 
gravity! 

The clearest way to illustrate this conclusion is to integrate away the metric ten- 
sor, i.e., to replace it in terms of its material sources, in order to construct the so- 
called "Fokker action". We give below a schematic derivation of its expression, 
taken from 1401 . Gauge-fixing subtleties are discussed notably in Appendix C of 
B3]| . We start from an action of the form 

5 = 5c(,[0]+5,natter[(T,^], (3) 

where denotes globally all fields participating in the gravitational interaction, and 
a denotes the matter sources. We also denote as ^[o] a solution of the field equation 
5S/8<P = for given sources a. Let us now define the Fokker action 

5Fokker[cT] = 5cf. [^[ff]] + ^matter [d, <i>[CT]] , (4) 

and show that it gives the correct equations of motion for matter (7. Indeed, its 
variational derivative reads 

\ J <P=<P[a] \ J <P=<P[a] 

where the second term of the right-hand side vanishes because [a] has been chosen 
as a solution of 5S/50 = 0. Therefore 5Spoi±er[(y]/S(y = does yield the correct 
equations of motion 5S[g, <P]/5g = for matter within the background = <l>[a] 
it consistently generates. The most important point to notice here is that the Fokker 
action (|4]l is not simply given by the matter action 5matter[<7, computed in the 
consistent background = 0[a]. Not only the a-dependence of this background 
must be taken into account when varying the Fokker action, but its definition (|4|i 
also depends crucially on the kinetic term (and the nonlinear dynamics) of the field, 

S4> [0[cy]\. 

To illustrate more vividly that the motion of massive bodies does depend on the 
dynamics of the gravitational field(s), let us give a diagrammatic representation of 
the above formal definition ^ of the Fokker action. We first introduce some sym- 
bols in Fig.[T] notably white blobs for matter sources and straight lines for field prop- 
agators. Using this notation, the original action ([3j may be translated as in Fig. |2] 
which actually defines the various vertices. In this figure, the numerical factors have 
been chosen to simplify the field equation satisfied by 0[(y], which takes the dia- 



a = O material sources 

□ <I> = • left hand side of the field equations 

Fig. 1 Diagrammatic repre- propagator of the fields (Green function) 

sentation of matter sources, ^ ^. 

fields, and their propagator. ^ ~ * ^ 
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kinetic term 
of tlie fields 



+ '^o[a] 



+ O 



fiee bodies linear interaction 
of matter and fields 



Yii: 

higher vertices 



Fig. 2 Diagrammatic representation of the full action {5} of the theory, expanded in powers of <P 
(black blobs). The first line corresponds to the field action S<p[<P], and the second one to the matter 
action 5matter[o'! (describing notably the matter-field interaction). 




Fig. 3 Diagrammatic representation of the field equation dS/S'P = satisfied by ^'[ct]. 



grammatic form of Fig. [3] This figure tells us how to replace any black blob (field <P) 
by a white blob (source a) plus higher corrections, in which one can again replace 
iteratively black blobs by white ones plus corrections. The Fokker action ^ is thus 
simply obtained by eliminating in such a way black blobs from Fig. El and the result 
is displayed in Fig.|4] This figure clearly shows that the dynamics of the field (i.e., 
the first line of Fig. |2]i does contribute to that of massive bodies. Indeed, if it had 
not been taken into account, the Newtonian interaction would have been twice too 



•SpotterW =So[o] + (i O O) + {j 

free bodies Newton ~~ 




-O 



O 



1PM 




2PM 



Fig. 4 Diagrammatic representation of the Fokker action j4j, which depends only on matter 
sources <J (white blobs). The dumbbell diagram labelled "Newton" represents the Newtonian inter- 
action <x G, together with all velocity-dependent relativistic corrections. The 3-blob diagrams la- 
belled "1PM" represent first post-Minkowskian corrections, i.e., the ff(G') post-Newtonian terms 
as well as their full velocity dependence. The 4-blob diagrams labelled "2PM" represent second 
post-Minkowskian corrections <x . 
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large, no "T" diagram would have appeared in Fig.H] and the numerical coefficients 
of all other diagrams would have also changed. 

For any theory of gravity, whose field dynamics is imposed by ^gravity, one may 
now explicitly compute the diagrams entering Fig. |4] Of course, any gauge invari- 
ance must be fixed in order to define the field propagator as the inverse of the 
quadratic kinetic term (first dumbbell diagram of Fig. |2]i. In the case of GR, one 
may for instance fix the harmonic gauge, and the first line of Fig. |4]translate^ as the 
well-known Einstein-Infeld-Hoffmann action describing the interaction of several 
massive bodies labelled A, B, . . . : 



Jpokker 



/ dtmAC^\ll-\\/c^ 

A •' 

1 \-i f J GniAmB 



1,9.9. 3 



E /.r^ l + 2^(vi + v^)-2^(v..v.) 

1 yPPN 

^ {nAB ■ va) (n^s • vb) + —5- (va - \b) 



1 [,,9^i!!^^2r''_,)^J}.), 

Here tab denotes the (instantaneous) distance between bodies A and B, Uab is the 
unit 3-vector pointing from B to A, va is the 3-velocity of body A, and a sum over 
B y^Ay^C allows B and C to be the same body. The first line of Eq. (|6]l, noted [cr] 
in Figs. |2] and m merely describes free bodies in special relativity. The second line 
of Eq. ^ describes the 2-body interaction, i.e., the dumbbell diagram of Fig.|4]that 
we labelled "Newton". Its lowest-order term is indeed the Newtonian gravitational 
potential, and we have also displayed its first post-Newtonian (IPN) corrections, 
of order ^(v^/c^). Finally, the last line of Eq. ^ corresponds to the "V" and "T" 
diagrams labelled "1PM" in Fig.lH computed here at their lowest (IPN) order. 

The two coefficients and y''™ entering Eq. ^ are simply equal to unity in 
GR. They were introduced by Eddington ifSOl to describe phenomenologically other 
possible theories of gravity, although he did not have any specific model in mind. It 
happens that the most natural alternatives to GR, scalar-tensor theories (that we will 
introduce in Sec.|4]below), do predict different values for these two parameters. This 
comes from the fact that massive bodies can exchange scalar particles in addition 
to the usual gravitons of GR. If we represent gravitons as curly lines and scalar 
fields as straight lines, the diagrams contributing to the Fokker action (|6]l are indeed 
displayed in Fig.|5] The four diagrams involving at least one scalar line contribute 
to change the values of the Eddington parameters and Y^™- We will give their 
explicit values in Eq. (|9]l below, but we refer to Ref. ll40l for their derivation from 
diagrammatic calculations. 

Besides j3^™ and y^™, many other parameters may actually be introduced to 
describe the most general behavior of massive bodies at the first post-Newtonian 



' One needs to compute the integrals represented by the various diagrams to derive expression (|6j. 
See Ref. 1401 for explicit diagrammatic calculations. 
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Fig. 5 Diagrams contributing 
to tiie A'-body action JSJ at 
IPN order, in scalar-tensor 
theories of gravity. Graviton 
and scalar exchanges are 
represented respectively as 
curly and straight lines. 



A B 



o- 



B 

-O 




order. Under reasonable assumptions, notably that the matter action takes the metric 
form ^matter [matter, g^v] and that the gravitational interaction does not involve any 
specific length scale, Nordtvedt and Will 11 141 showed that 8 extra parameters are a 
priori possible, in addition to Eddington's jS^™ and y''™. However, these 8 param- 
eters vanish both in GR and in scalar-tensor gravity, therefore we will not introduce 
them in these lecture notes. A detailed presentation is available in the book 111 121 . 



4 Scalar-tensor theories of gravity 

Among alternative theories of gravity, those which involve scalar partners to the 
graviton are privileged for several reasons. Not only their existence is predicted in all 
extra-dimensional theories, but they also play a crucial role in modern cosmology (in 
particular during the accelerated expansion phases of the Universe). They are above 
all consistent field theories, with a well-posed Cauchy problem, and they respect 
most of GR's symmetries (notably conservation laws and local Lorentz invariance 
even if a subsystem is influenced by external masses). To simplify the discussion, 
we will focus on models involving a single scalar field, although the study of tensor- 
multi-scalar theories can also be done in great detail l38l . We will thus consider the 
class of theories defined by the action lfT9l [85l II 101 

5= J d'^xy/^ {R* ^Igi^^d^cpdycp) +5matter [matter, g = {(p) g*^^] . 

(7) 

A potential V{(p) may also be considered in this action, and is actually crucial in 
cosmology, but we will study here solar-system-size effects and assume that the 
scalar-field mass (and other self-interactions described by V{(p)) is small enough to 
be negligible at this scale. The physical metric gf^y, to which matter is universally 
coupled (and which defines thus the lengths and times measured by material rods 
and clocks), is the product of the Einstein metric y (whose kinetic term is the 
Einstein-Hilbert action) and a function A^((p) characterizing how matter is coupled 
to the scalar field. It will be convenient to expand it around the background value (po 
of the scalar field far from any massive body, as 



lnA((p) =lnA((j!)o) + 0!o(9-<Po) + 2i3o(<p-<Po)" + ^(<P-<Po)^ (8) 
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where oCq defines the Hnear couphng constant of matter to scalar excitations, j3o its 
quadratic couphng to two scalar lines, etc. 



4.1 Weak-field predictions 

Newtonian and post-Newtonian predictions depend only on these first two cou- 
pling constants, Oq and j3o. For instance, the effective gravitational constant be- 
tween two bodies is not given by the bare constant G* entering action (|7]i, but by 
G = G* (1 + ttg), in which a contribution G* comes from the exchange of a (spin-2) 
graviton whereas G* Oq is due to the exchange of a (spin-0) scalar field, each matter- 
scalar vertex bringing a factor Oq. The first line of Fig. |5] gives a diagrammatic il- 
lustration of this sum. [Actually, the value of a gravitational constant depends on 
the chosen units, and the expression G = G*(l + aij) corresponds to the "Einstein- 
frame" representation used to write action (|7]i. An extra factor Aq = A(^o)^ enters 
when using the physical metric g^^ ~ {^)g]iv 1° define observable quantities, and 
the actual gravitational constant which is measured reads G*Aq(1 + Ctg ). No such ex- 
tra factors Aq enter the computation of dimensionless observable quantities, like the 
Eddington parameters jS''™ and y''™.] Two kinds of post-Newtonian corrections 
enter the Fokker action (|6]i: velocity-dependent terms in the 2-body interaction (first 
line of Fig. Is), which involve the parameter 7''™, and the lowest-order 3-body inter- 
actions (second line of Fig.|5]l, which involve j3''™. Diagrammatic calculations ll42l 
or more standard techniques 11 121 [38l can be used to compute their expressions in 
scalar- tensor theories: 

,PPN _ 1 2*^0 oPPN _ , , 1 "0ft"0 /Qx 



^ Lunar Laser Ranging 



Fig. 6 Solar-system con- 
straints on the post-Newtonian 
parameters jS*"™ and "f™. 
The allowed region is the tiny 
intersection of the "Lunar 
Laser Ranging" strip with 
the horizontal bold line la- 
belled "Cassini". General 
relativity, corresponding to 

^PPN ^ ^PN ^ 1^ 

tent with all tests. 



1.004 



1.002 




0.998 .r-W-; X- 



0.996 
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Fig. 7 Solar-system con- 
straints on the matter- 
scalar coupling function 
lnA((p)/Ao = ao{(p-(po) + 
5A)(<P-<Po)^ + -- - The al- 
lowed region is the dark 
grey horizontal strip. The 
vertical axis (jSq = 0) corre- 
sponds to Brans-Dicke theory 
dUlTOlIll] with a param- 
eter 2(Bbd + 3 = 1/al The 
horizontal axis (Oq = 0) cor- 
responds to theories which are 
perturbatively equivalent to 
GR, i.e., which predict strictly 
no deviation from it (at any 
order 1/c") in the weak- field 
conditions of the solar system. 



Here again, the factor comes from the exchange of a scalar particle between two 
bodies, whereas OojSoCto comes from a scalar exchange between three bodies (cf. the 
purely scalar "V" diagram of Fig.|5]l. 

Several solar-system observations tightly constrain these post-Newtonian param- 
eters to be close to 1, i.e., their general relativistic values. The main ones are Mer- 
cury's perihelion advance 195)1 , Lunar Laser Ranging (which allows us to test the 
so-called Nordtvedt effect, i.e., whether there is a difference between the Earth's 
and the Moon's accelerations towards the Sun) 111 151 . and experiments involving 
the propagation of light in the curved spacetime of the solar system (by order of 
increasing accuracy: radar echo delay between the Earth and Mars, light deflection 
measured by Very Long Baseline Interferometry over the whole celestial sphere 
ll96ll . and time-delay variation to the Cassini spacecraft near solar conjunction ||20| ). 
These post-Newtonian constraints are summarized in Fig. |6] and the conclusion is 
that GR is basically the only theory consistent with weak-field experiments. How- 
ever, when translated in terms of the linear and quadratic coupling constants Oq and 
/3o of matter to the scalar field, the same solar-system constraints take the shape of 
Fig. Ill Therefore, the linear coupling constant |ao| must be smaller than 3 x 10^^, 
but we do not have any significant constraint on /3o [nor any higher-order vertex 
entering expansion dHJ]. 



|aol O 




4.2 Strong-field predictions 

A qualitatively different class of constraints is obtained by studying scalar-tensor 
theories in the strong-field regime, i.e., near compact bodies whose radius R is not 
extremely large with respect to their Schwarzschild radius 2Gm/c^. This is notably 
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the case when considering neutron stars, whose ratio Gm/Rc^ ^ 0.2 is not far from 
the theoretical maximum of 0.5 for black holes. Although the metric is very signifi- 
cantly different from the flat one inside such compact bodies and in their immediate 
vicinity, their orbital velocity in a binary system may nevertheless be small enough 
to perform a consistent expansion in powers of v^/c^ ~ Gm/rc^, where r denotes 
the interbody distance (as opposed to their radius R). This is usually called a post- 
Keplerian expansiorH In such a case, one can show 11 121 [38l that the predictions 
of scalar-tensor theories are similar to those of weak-field conditions, with the only 
difference that the matter-scalar coupling constant Oq and jSo are replaced by body- 
dependent quantities, say Ua and jS^ for a body labelled A. For instance, the effective 
gravitational constant describing the lowest-order attraction between two compact 
bodies A and B reads now Gab = Gh,(1 + aAOLs), instead of the weak-field expres- 
sion G = G*(l + Oq) mentioned above. Similarly, the post-Newtonian parameters 
/3^™ and 7''™ are replaced by body-dependent ones /3g^ and Jab, taking the same 
forms as in Eq. (|9]l but where Oq is replaced by UaOLb, and (Xo^oOk) by Ub^aOIc (see 
Ref. 1381 for precise expressions). All post-Keplerian effects can thus be derived 
straightforwardly, in a similar way as in the solar system. In addition to these pre- 
dictions, one may also compute the energy loss due to the emission of gravitational 
waves by a binary system. It takes the schematic forrr0 



Energy ^ f Quadrupole ^ / j_ 

flux "\ c5 yc'' 



spin 2 



I Monopole ^ Jl ^^ ^ ^ Digole ^ Quadrupole + ^ i_ | , (iq) 



spin 

where the first line comes from the emission of usual (spin-2) gravitons, and the 
second one from the emission of scalar (spin-0) waves. Note that the dipolar term is 
of order 0'{l/r'), generically much larger than the standard 6'[\/r') quadrupole of 
GR. As expected for a dipole, it vanishes when considering a perfectly symmetrical 
binary system, because there is no longer any privileged spatial orientation. Its pre- 
cise calculation 11 121 [38l shows indeed that it involves the difference of the scalar 
charges of the two bodies, and is actually proportional to (a^ — OLb)^- Although the 
monopolar term is a priori of the even larger order ^( 1 /c) for bodies which are not 
at equilibrium (e.g. collapsing or exploding stars), it reduces to order ff{llc^) for 
usual bodies, as displayed in Eq. ( fTol i. because of the conservation of their scalar 
charge. 



^ Two different (though related) meanings of "post-Keplerian" exist in the literature. We are here 
considering a post-Keplerian expansion in powers of v^|.|^jj^|/c"^, while keeping the full nonper- 
turbative dependence in the gravitational self-energy Gm/Rc^ . On the other hand, post-Keplerian 
deviations mean relativistic effects modifying the lowest-order Keplerian motion, like those de- 
scribed in Sec. |4.3| below. Only this latter meaning is used in GR, because its strong equivalence 
principle implies that the internal structure of a body does not influence its motion up to order 
ff{l/c^^), as recalled in Sec.|5]below. 

^ The precise definitions of these multipoles and their explicit expressions may be found for in- 
stance in Sec. 6 of Ref. 1381. 
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Fig. 8 Scalar charge (Za 
of a neutron star versus its 
baryonic mass mA, for the 
model A{(p) = exp(-3(p^), 
i.e., jSo = —6. The solid line 
corresponds to a small value 
of Oo (namely, the VLBI 
bound of Fig. |7), and the 
dashed line to ao = 0. The 
dotted lines correspond to 
unstable configurations of the 
star. 
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The only remaining difficulty, to derive the predictions of scalar-tensor theories 
in the strong-field regime, is to compute the body-dependent coupling constants Ua 
and This can be done thanks to numerical integrations of the field equations in- 
side the bodies, as explained in ||39ll4Tll42l . For negative values of the parameter j3o 
entering expansion ([8]l, one shows, both analytically and numerically, that nonper- 
turbative effects occur beyond a critical compactness Gm/Rc^ depending on Pq. For 
instance, for jSo ~ —6, the linear coupling constant of a neutron star to the scalar 
field takes the values displayed in Fig. [8] One sees that even if oco is vanishingly 
small in the background (and thereby in the solar system), neutron stars can develop 
an order-one coupling constant to the scalar field. Their physics and their orbital 
motion can thus differ significantly from the predictions of GR, although the scalar 
field may have strictly no effect in the solar system. 



4.3 Binary-pulsar tests 

Now that we know how to compute the predictions of scalar-tensor theories even 
in strong-field conditions, how may we test them? It happens that nature has pro- 
vided us with fantastic objects called pulsars. These are neutron stars (thereby very 
compact objects, Gm/Rc" ~ 0.2) which are rapidly rotating and highly magnetized, 
and which emit a beam of radio waves like lighthouses. They can thus be consid- 
ered as natural clocks, and the oldest pulsars are indeed very stable ones. Therefore, 
a pulsar A orbiting a companion B is a moving clock, the best tool that one could 
dream of to test a relativistic theory. Indeed, by precisely timing its pulse arrivals, 
one gets a stroboscopic information on its orbit, and one can measure several rel- 
ativistic effects. Such effects do depend on the two masses m^, ms, which are not 
directly measurable. However, two different effects suffice to determine them, and 
a third relativistic observable then gives a test of the theory. 

For instance, in the case of the famous Hulse-Taylor binary pulsar PSRB1913H-16 
111 1 II . three relativistic parameters have been determined with great accuracy: (i) the 
Einstein time delay parameter Yt, which combines the second-order Doppler ef- 



Title Suppressed Due to Excessive Length 



13 



feet (0= v^/2c^) together with the redshift due to the companion (0= GniQ/r/^QC^); 
(ii) the periastron advance 6) v^/c^); and (iii) the rate of change of the orbital 
period, P, caused by gravitational radiation damping (0= v^/c^ in GR, but of order 
v^/c^ in scalar-tensor theories; see Eq. (fTOli). The same parameters have also been 
measured for the neutron star-white dwarf binary PSR Jl 141— 6545, but with less 
accuracy ifTOl fTTI . In addition to these three parameters, (iv) the "range" (global 
factor GrriB /r') and (v) "shape" (time dependence) of the Shapiro time delay have 
also been determined for two other binary pulsars, PSR B1534H-12 HI 021 and PSR 
J0737— 3039 Il30ll75ll72l . The latter system is particularly interesting because both 
bodies have been detected as pulsars. Since their independent timing gives us the 
(projected) size of their respective orbits, the ratio of these sizes provides a direct 
measure of (vi) the mass ratio mp^/m-^ « 1.07. In other words, 6 relativistic param- 
eters have been measured for the double pulsar PSR J0737— 3039. After using two 
of them to determine the masses and me, this system thereby provides 6 — 2=4 
tests of relativistic gravity in strong-field conditions. 

The clearest way to illustrate these tests is to plot the various experimental con- 
straints in the mass plane [mx^m-^), for a given theory of gravity. Any theory indeed 
predicts the expressions of the various timing parameters in terms of these unknown 
masses and other Keplerian observables, such as the orbital period and the eccentric- 
ity. The equations predictions {m/^,ms) = observed values thereby define different 
curves in the mass plane, or rather different strips if one takes into account exper- 
imental errors. If these strips have a common intersection, there exists a pair of 
masses which is consistent with all observables, and the theory is confirmed. On the 
other hand, if the strips do not meet simultaneously, the theory is ruled out. Figure|9] 
displays this mass plane for the Hulse-Taylor binary pulsar. Its left panel shows 




0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 



Fig. 9 Mass plane (iha = pulsar, = companion) of the Hulse-Taylor binary pulsar PSR 
B 19 13+16 in general relativity (left panel) and for a scalar-tensor theory with /3o = — 6 (right 
panel). The widths of the lines are larger than 1(7 error bars. While GR passes the test with flying 
colors, the value j8o = —6 is ruled out. [Actually, the right panel is plotted for a specific small value 
of |c<o| Ri 10^^, but the three curves keep similar shapes whatever Oq, even vanishingly small, and 
they never have any common intersection for j8o ~ —6.] 
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Fig. 10 Same theory plane 
(|ao|,A)) as in Fig. El but 
taking now into account the 
constraints imposed by the 
Hulse-Taylor binary pulsar 
(bold line). LLR stands as 
before for "lunar laser rang- 
ing" and VLBI for "very long 
baseline interferometry". The 
allowed region is the dark 
grey one. While solar system 
tests impose a small value 
of |c<o|, binary pulsars im- 
pose the orthogonal constraint 
j8o > -4.5. 
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that GR is superbly consistent with these data. Its right panel illustrates that the 
three strips can be significantly deformed in scalar-tensor theories, because scalar 
exchanges between the pulsar and its companion modify all theoretical predictions. 
In the displayed case, corresponding to a quadratic matter-scalar coupling constant 
j3o = —6 (as in Fig.[8]i, the strips do not meet simultaneously and the theory is thus 
excluded. On the contrary, they may have a common intersection in other scalar- 
tensor theories, even if it does not correspond to the same values of the masses 
and nis that were consistent with GR. The allowed region of the theory space 
(|ao|,j3o) is displayed inFig.fTOl 

As mentioned above, several other relativistic binary pulsars are presently known, 
and Fig. [TT] displays their simultaneous constraints on this theory plane BT[|42ll54l 
[36l|43|. To clarify this plot, we have used a logarithmic scale for the vertical (|ao|) 
axis. The drawback is that GR, corresponding to Oo = j3o = 0, is sent down to in- 
finity, but the important point to recall is that it does lie within the allowed (grey) 
region. Figures[TO]and[TT]illustrate vividly the qualitative difference between solar- 
system and binary -pulsar observations. Indeed, while weak-field tests only constrain 
the linear coupling constant | oco | to be small without giving much information about 
the quadratic coupling constant )3o, binary pulsars impose /3o > —4.5 even for a van- 
ishingly small Oq. This constraint is due to the spontaneous scalarization of neutron 
stars which occurs when — j3o is large enough, and which was illustrated in Fig. [8] 
above. Equations (|9]l allow us to rewrite this inequality in terms of the Eddington 
parameters and y^™, which are both consistent with 1 in the solar system. 
One finds 

gPPN _ 1 

^p™^<1.1. (11) 

The singular (0/0) nature of this ratio underlines why such a conclusion could not 
be obtained in weak-field experiments. 
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Fig. 11 Same tlieory plane 
as in Fig. [TO] but now with 
a logarithmic scale for the 
linear matter-scalar coupling 
constant | Oo | . This plot dis- 
plays solar-system and several 
binary-pulsar constraints, us- 
ing the latest published data. 
The allowed region is shaded. 
The curve labeled SEP corre- 
sponds to tests of the "strong 
equivalence principle" using 
a set of neutron star-white 
dwarf low-eccentricity bi- 
naries. The dashed lines 
corresponds to expected con- 
straints if we find a relativistic 
pulsar-black hole (PSR-BH) 
binary, or when gravity-wave 
antennas detect the coales- 
cence of double-neutron star 
(NS-NS) or neutron star-black 
hole (NS-BH) binaries. 
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4.4 Black holes in scalar-tensor gravity 

In Fig. [TT]is displayed the order of magnitude of the tight constraint that can be 
expected if we are lucky enough to detect a relativistic pulsar-black hole binary. 
However, let us comment on the behavior of black holes in scalar-tensor theories. 
Since we saw in Fig. [8] above that nonperturbative effects can occur in strong-field 
conditions, one should a priori expect even larger deviations from GR for black 
holes (extreme compactness Gm/Rc^ = 0.5) than for neutron stars (large compact- 
ness Gm/Rc^ - 0.2). However the so-called no-hair theorem lfT3ll62ll27ll60ll59ll38l 
shows that black holes must have a strictly vanishing scalar charge, Ubh = 0. The 
basic idea is that otherwise the scalar field <p would diverge at the horizon, and 
this would be an unphysical solution. A first consequence is that a collapsing star 
must radiate away its scalar charge when forming a black hole. This is related to 
the generically large ff{\/c) monopolar radiation of scalar waves predicted for non- 
equilibrium configurations, as discussed below Eq. ( fTOl ). But the second crucial con- 
sequence is that black holes, once formed and stabilized, are not coupled at all to the 
scalar field, and therefore behave exactly as in GR: They generate the same solution 
for the Einstein metric g^y, do not excite the scalar field <p, and move within the 
curved geometry of g*^^ as in GR. The conclusion is therefore that there is strictly 
no observable scalar-field effect in a binary black-hole system. 

Of course, there do exist significant perturbations caused by the scalar field dur- 
ing the short time of the black-hole formation or when it captures a star (see e.g. 
II3TI ). because of the emission of a generically large amount of energy via scalar 
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waves. Similarly, if one assumes that there exists a non-constant background of 
scalar field (po(x), then its own energy momentum tensor Tfiv contributes to the cur- 
vature of the Einstein metric g^y, and even black holes would thus indirectly feel its 
presence. However, this would not be a consequence of the modification of gravity 
itself, but of the assumption of a non-trivial background T^v Even within pure GR, 
one could also have assumed that black holes move within a non-trivial background, 
caused for instance by the presence of dark matter or some large gravitational waves, 
and one would predict then a different motion than in a trivial background. One 
should thus qualify our conclusion above: Given some precise boundary conditions, 
for instance asymptotic flatness and no incoming radiation, black holes at equilib- 
rium behave exactly as in GR. 

It remains now to understand how a pulsar-black hole binary can allow us to 
constrain scalar-tensor theories as in Fig. [TT] The reason is that one of the two 
bodies is not a black hole, but a neutron star with scalar charge Ua ^ 0. Let us 
even recall that massive enough neutron stars can develop order-one scalar charges 
as in Fig. [8] even if matter did not feel at all the scalar field in the solar system 
(Oo = 0). A pulsar-black hole binary must therefore emit a large amount of dipolar 
waves o= (ttA — O-bh)^ /c^ = ^a/^^ — ^(l/*^^)' given by Eq. ( fTol i. and this can 
be several orders of magnitude larger than the usual &{\/c^) quadrupolar radiation 
predicted by GR. Any pulsar-black hole binary whose variation of the orbital period, 
P, is consistent with GR will thus tightly constrain scalar-tensor models. 



5 Extended bodies 

In this section, we will discuss finite-size effects on the motion of massive bodies, 
in GR and in scalar-tensor theories ll86l l42l 1611 . It is still convenient to describe the 
position of such extended bodies by using one point in their interior, for instance 
their approximate center of mass. In other words, we will skeletonize the extended 
body's worldtube as a unique worldline, say x^. However, the action (|2]i describing 
the motion of a point particle cannot remain valid to all orders, because the metric 
g^v(jcieft) on one side of the body is not strictly the same as g/jv (-bright) on the other 
side. By expanding the metric around its value at Xc, it is thus clear that an effective 
action describing the motion of an extended body must depend on derivatives of 
the metric, and since we wish to construct a covariant expression, such derivatives 
must be built from contractions of the curvature tensor, its covariant derivatives, 
or any product of them. Moreover, if the body is nonspinning and spherical when 
isolated, it does not have any privileged direction, and the only 4-vector available to 
contract possible free indices is the 4-velocity = dx^ jds of the point Xc- The first 
couplings to curvature that one may think of are thus of the form IMI 

■^extended body — 5point particle + / (^1^ + kiRy^yUy-u' ^ ) Cc/s, (12) 
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where and k2 denote body-dependent form factors (no term Ri^ypaii^u^u^u" is 
written since it vanishes because of the symmetry properties of the Riemann ten- 
sor). However, as recalled in 14211241 . any perturbative contribution to an action, say 
5[v/], which is proportional to its lowest-order field equations, 5S/5\l/, is equivalent 
to a local field redefinition. Indeed, S[\(/ + £] = S[\(f] + e5S/5\if+ ff{e^), where the 
small quantity e may be itself a functional of the field Therefore, up to local 
redefinitions of the worldline Xc and the metric inside the body, which do not 
change the observable effects encoded in the metric gfivixoha) at the observer's lo- 
cation, one may replace R and /J^v in action ( fT2] i above by their sources. Outside the 
extended body, such couplings to R or R/jy have thus strictly no observable effect 
(see ll24l for a discussion of the field redefinitions inside the body even when it is 
skeletonized as a point particle). 

The first finite-size observable effects must therefore involve the Weyl tensor 
Cfivpo, which does not vanish outside the body. Because of the symmetry properties 
of this tensor, the lowest-order terms one may construct from it take the form lIMl 



By comparing the dimensions of S'point particle — ~S mcds and these couplings to 
the Weyl tensor, we expect kC^ ^ m and therefore k ^ mR'j^ (where k denotes any 
of the form factors k^,, ^4 or ^5, and Re means the radius of the Extended body, 
here written with an index to avoid a confusion with the curvature scalar). In con- 
clusion, for a compact body such that Gm/Rec^ ~ 1, we expect the form factors 
k to be of order m{Gm/c^Y, and the corrections to its motion will be proportional 
to kc^C'^ = kc^i&{l /c^) = /c'"). Therefore, finite size effects start at the fifth 
post-Newtonian (5PN) order, consistently with the Newtonian reasoning. Indeed, if 
one considers two bodies A and B which are spherical when isolated, one can show 
that B is deformed as an ellipsoid by the tidal forces caused by A, and this defor- 
mation induces an extra force ^ [GniA/ r\g){RB / rAs)^ felt by A. If B is a compact 
body, i.e., GniB /Rrc^ ^ 1, one concludes that the extra force felt by body A is of or- 
der /c^"). We recovered above the same conclusion within GR thanks to a very 
simple dimensional argument. Note however that the explicit calculation is more 
involved than our estimate of the order of magnitude, because one must take into 
account the Weyl tensor generated by the two (or more) bodies of the system. 

Let us now follow the same dimensional reasoning within scalar-tensor theories 
of gravity. We assumed in Eq. (|7]i that matter is coupled to the physical metric gfxv = 
A^((p)g^y, therefore the action of a point particle reads 



5point particle = - / mc^ -g^ydx^dx'^ = - / A{^)mc J -g*^^dx^ dx'^ . (14) 



When studying the motion of such a point particle in the Einstein metric g*^,, it 
behaves thus as if it had a scalar-dependent mass m*((p) = A{<p)m. [Incidentally, 
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this provides us with another definition of the scalar charge, a = t/lnm* {(p) /d(p.] If 
we consider now an extended body in scalar-tensor gravity, the function (9) must 
be replaced by afunctional, i.e., it can depend on the (multiple) derivatives of both 
the metric and the scalar field. Let us restrict again our discussion to nonspinning 
bodies which are spherical when isolated. Then the derivative expansion of this mass 
functional can be written as ll42ll 

mW,8liv] = m{(p)+I{(p)R*+J{(p)Rl^u^ul+K{(p)n*(p + L{(p)Vl,dy(pu'iul 

+M{(p)d^(pdv(pu^ul +N{(p)g^''d^(pdv(p + ■■■, (15) 

where star-indices mean that we are using the Einstein metric g*^y to compute cur- 
vature tensors, covariant derivatives and the unit velocity = dx^^ /ds*, and where 
/((p), . . . , N{(p) are body-dependent form factors encoding how the extended body 
feels the various second derivatives of the fields (p and g*^^. Fortunately, by using 
as in GR the lowest-order field equations together with local field and worldline 
redefinitions, this long expression reduces to a mere 

m[(p,g*^y] = m{(p) +A^new(^)^*^^|i<P^v<P + higherpost-Keplerian orders, (16) 

where A^new(<p) is a linear combination of the previous N{(p), L{(p) and I{(p). We 
already notice that the lowest-order finite-size effects involve only two derivatives 
of the scalar field, whereas action (T3[ in GR involved the square of a second deriva- 
tive of the metric. As above, let us invoke dimensional analysis to deduce that 
A^new ~ mRj^, where Re still denotes the radius of the extended body. This is consis- 
tent with the actual calculation performed in ll86l for weakly self-gravitating bodies, 
where it was proven that A^new = 5)^0 x (Inertia moment). For compact (therefore 
strongly self-gravitating) bodies, such that Gm/Rec^ ^ 1, one thus expects observ- 
able effects of order Nc-{d^f = ff{\/c^)R\ = (^(l/c*^), i.e., of the third post- 
Newtonian (3PN) order Finite-size effects are thus much larger in scalar-tensor the- 
oriefl than in GR, where they were of the 5PN order, (^(l/c'''). Moreover, when 
nonperturbative strong-field effects develop as in Fig.|8]above, one actually expects 
that A^new ~ X (Inertia moment) can be of order unity, because the coupling con- 
stant Pe of the body to two scalar lines can be extremely large with respect to the 
bare jSo. Therefore, finite-size effects in scalar-tensor gravity should actually be con- 
sidered as first post-Keplerian, ^(^'orbitai/'-^)' perturbations of the motion, as com- 
pared to the fifth post-Newtonian order in GR. 



These larger finite-size effects are due to the fact that a spin-0 scalar field can couple to the 
spherical inertia moment of a body, contrary to a spin-2 graviton. They should not be confused 
with the violation of the strong equivalence principle, which also occurs in scalar-tensor theories 
because all form factors m{(p), N{(p), . . . depend on the body's self-energy. Regardless of its finite 
size, the motion of a self-gravitating body in a unifonn exterior gravitational field depends thus on 
its internal structure. 
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6 Modified Newtonian dynamics 

The existence of dark matter (a pressureless and non-interacting fluid detected only 
by its gravitational influence) is suggested by several observations. For instance, 
type-la supernova data Ill041l7n[89l l7l are consistent with a present acceleration of 
the expansion of the Universe, and tell us that the dark energy density should be of 
order Qj^ k, 0.7. On the other hand, the position of the first acoustic peak of the cos- 
mic microwave background spectrum 11001 IIOII is consistent with a spatially flat 
Universe, i.e., Q.j^ +i2matter ~ 1- Combining these two pieces of information, one 
thus deduce that the matter density should be f2matter ~ 0.3, a value at least one order 
of magnitude larger than all our estimates of baryonic matter in the Universe (for 
instance X2baryons ~ 0.04 derived from Big-Bang nucleosynthesis). Therefore, most 
of the cosmological matter should be non-baryonic, i.e., "dark" because it does not 
interact significantly with photons. Another, independent, evidence for the existence 
of dark matter is the flat rotation curves of clusters and galaxies [90| : The velocities 
of outer stars tend toward a constant value (depending on the galaxy or cluster), in- 
stead of going asymptotically to zero as expected in Newtonian theory (recall that 
Neptune is much slower than Mercury in the solar system). If Newton's law is as- 
sumed to be valid, such nonvanishing asymptotic velocities imply the existence of 
much more matter than within the stars and the gas. Independently of these experi- 
mental evidences, we also have many theoretical candidates for dark matter, notably 
the class of neutralinos occurring in supersymmetric theories (see e.g. fZSl ). and 
numerical simulations of structure formation have obtained great successes while 
incorporating dark matter (see e.g. 163)1 ). 

However, this unknown fluid might actually be an artifact of our interpretation of 
experimental data with a Newtonian viewpoint. It is thus worth examining whether 
the gravitational l/r^ law could be modified at large distances, instead of invoking 
the existence of dark matter. In 1983, Milgrom realized that galaxy rotation curves 
could be fitted with a very simple recipe, that he called Modified Newtonian Dy- 
namics (MOND) 1781 . It does not involve any mass scale nor distance scale, but an 
acceleration scale denoted as aq [not to be confused with the linear matter-scalar 
coupling constant Oq, defined in Eq. (H)) above]. Milgrom assumed that the acceler- 
ation fl of a test particle caused by a mass M should read 

GM 

a— aN — — ita>ao, 

(17) 

^JGMao ^ 

a ~ y/aoajM = it a < ao, 

where afj denotes the usual Newtonian expression. This phenomenological law hap- 
pens to fit remarkably well galaxy rotation curves |j94l . for a universal constant 
ao ~ 1.2 X lO^'^m.s^^. [However, galaxy clusters require either another value of 
this constant, or some amount of dark matter, for instance in the form of massive 
neutrinos.] Moreover, it automatically recovers the Tully-Fisher law 11051 vt, °^ M, 
where M denotes the baryonic mass of a galaxy, and Voo the asymptotic velocity of 
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visible matter in its outer region. The MOND assumption ( fTTj i would also explain 
in an obvious way why dark matter profiles seem to be tightly correlated to the 
baryonic ones 1761 . 

However, reproducing the simple law ( fTTb in a consistent relativistic field the- 
ory happens to be quite difficult. As mentioned in Sec. |2] above, Milgrom explored 
"modified inertia" models 1791 [80l . in which the action of a point particle is as- 
sumed to depend nonlocally on all the time derivatives d"x/dt" of its position x. 
In these lecture notes, we focus on local field theories which "modify gravity", i.e., 
which assume that the kinetic term of the metric g^v (to which matter is universally 
coupled) is not the standard Einstein-Hilbert action ([T]i. 



6.1 Mass -dependent models? 

It is actually very easy to devise a model which predicts a force 0= 1 /r. Indeed, let us 
consider a mere scalar field (p in flat spacetime, with a potential V{(p) = —la^e^'"^, 
where a and b are two constants. In a static and spherically symmetric situation, its 
field equation A(p ^ ^'i^P) then gives the obvious solution (p — {2/b)ln{abr). If we 
now assume that matter is linearly coupled to (p, it will feel a force di(p, i.e., the 
1/r law we are looking for However, this simple model presents two very serious 
problems. The first one is that the potential V{(p) ~ —la^e^^f is unbounded from 
below, and therefore that the theory is unstable. Actually, we will see that stability is 
indeed a generic difficulty of all models trying to reproduce the MOND dynamics. 
The second problem of this naive model is that it predicts a constant coefficient 2//? 
for the solution of the scalar field, instead of the factor \/M entering the second line 
of ([TtI i. It happens that many models proposed in the literature do behave in the 
same way, although their more complicated writings hide the problerr0. The trick 
used by the corresponding authors is merely to set 1 / VM in the action of the 
theory. In other words, they are considering a different theory for each galaxy Ml It 
should be noted that the mass of an object is not a local quantity: It is the integral of 
the matter density over a particular region, whose symmetry plays also an important 
role. One of the most difficult steps in building a consistent theory of MOND is thus 
precisely to be able to predict this factor VM- Moreover, if one defined a potential 
V{(p) as above, in terms of an integral of the matter density giving access to M, then 
it would mean that matter is coupled to ^ in a highly nonlinear and nonlocal way, 
therefore the force it would feel would be much more complicated than the naive 
gradient di(p assumed above. In the following, we will consider local field theories 
whose actions do not depend on M, but only on the constants G, c and oq. 



^ See Sec. II.B of Ref. [291 for a critical discussion of various sucli mass-dependent models. 
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6.2 Aquadratic Lagrangians or k-essence 

One of the most promising frameworks to reproduce the MOND dynamics is a gen- 
eralization of the scalar- tensor theories we considered in Sec. |4] above. Their action 
takes the form |fT4l[T8ll92l[T5l[T6ll93l 

S = / d'x^ {R* ~2f{g^'d,(pd,(p) } 

+5matter [matter;g^v = A-{(p)gl,y+B{(p)Uf,Uv] . (18) 

The first crucial difference with action (|7J is that the kinetic term of the scalar field 
is now a function of the standard quadratic term {df^cp)^. Reference ifTTl showed 
that this suffices to reproduce the MOND law ( [TtI i. including the important \/M 
coefficient, provided /(x) 0= x^/'^ for small values of x (MOND regime) and /(x) 
X for large x (Newtonian regime). Such aquadratic kinetic terms have also been 
analyzed later in the cosmological context, under the names of k-inflation Q or 
k-essence Il32l l6l (the letter k meaning that their dynamics is kinetic dominated). 

As in action matter is assumed to be coupled to the scalar field via the func- 
tion A^{(p) entering the definition of the physical metric gi^y. This ensures that test 
particles will undergo an extra acceleration caused by the scalar field. But the second 
difference with (|7]i is the presence of the non-conformal term B{<p)UfiUv entering 
the definition of Here [/^ is a vector field, which can have either its own kinetic 
term, or can be simply chosen as C/^ — dp^cp. This term is necessary to reproduce 
the light deflection caused by galaxies or clusters, which happens to be consistent 
with the prediction of GR in presence of a dark matter halo ll56l 1771 [121 . Indeed, 
light is totally insensitive to the conformal factor A^((p) relating the physical and 
Einstein metrics in action (|7]i- The simplest way to understand it is to recall that a 
photon's null geodesic satisfies ds^ = A~{(p)dsl = dsl ~ 0, therefore the 
photon propagates in the Einstein metric y without feeling the presence of the 
scalar field. [To be more precise, the photon does feel indirectly the presence of the 
scalar field via the extra curvature of g*^y caused by its energy-momentum tensor, 
but this is a higher post-Newtonian effect.] The only way to impose that the MOND 
potential (encoded in the scalar field (p) deflect Ught is thus to relate the physical 
metric to the Einstein one in a non-conformal way, as in ( fTSI ). This idea dates 
back to Ni's "stratified" theory of gravity ll83lll 121 . 

Obviously, several conditions must be imposed on the functions entering action 
(fTSb to warrant that the theory is stable and that it has a well-posed Cauchy problem. 
For instance, it is clear that /(x) = x defines a standard kinetic term for the scalar 
field, whereas /(x) = — x would define a negative-energy (ghost) mode. In order 
for the Hamiltonian to be bounded by below and for the scalar-field equations to be 
hyperbolic, one can actually show that the function / must satisfy the two conditions 

Eiiaii 

Vx, /'(x)>0, and Vx, 2x/"(x) +/'(x) > 0. (19) 
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One also notices that gravitons are faster than scalar particles when f"{x) < 0, and 
slower when f"{x) > 0. It has been argued in Ul that gravitons should be the fastest 
modes, otherwise causal paradoxes can be constructed. However, Refs. ||28l [8l l29l 
concluded that conditions ( fT9] l suffice for causality to be preserved, because they 
ensure that the widest causal cone always remains a cone, and never opens totally. 

The analogues of conditions ( fT9] l become much more complicated within matter, 
where the two other functions A{<p) and B{<p) of action (fTsl l also enter the game. 
Moreover, one also needs to ensure that the matter field equations always remain 
hyperbolic. We refer to ||29l for a discussion of these issues. 



6.3 Difficulties 

Although the class of relativistic models ( fTSl l is the most promising one to repro- 
duce the MOND phenomenology, it has anyway a long list of difficulties. Some of 
them can be solved, but at the price of complicated and unnatural Lagrangians. For 
instance, we mentioned above the problem of light deflection, which is too small in 
conformally coupled scalar-tensor theories (|7]i, therefore one needed to introduce a 
vector field U^i in action (fTsT l. One may also notice that such models are not very 
predictive, since it would have been possible to predict fully different lensing and 
rotation curves. One may thus consider them as fine-tuned fits rather than fundamen- 
tal theories imposed by some deep symmetry principles. The most famous model, 
called TeVeS (for Tensor- Vector-Scalar) llT5l[T6l , also presents some discontinuous 
functions, and does not allow to pass smoothly from a time evolution (cosmology) 
to a spatial dependence (local physics in the vicinity of a galaxy). However, some 
cures are possible ||29l , although they again involve rather unnatural refinements. 

One peculiarity of the TeVeS model is that its author imposed that gravitons and 
scalar particles are slower than photons (a priori to avoid causal paradoxes due 
to superluminal propagation, assuming that light is a more fundamental field than 
gravity). But Refs. 1821 l53l proved that in such a case, high-energy cosmic rays 
would rapidly lose their energy by Cherenkov radiation of gravitational waves, and 
this would be inconsistent with their observation on Earth. However, a very simple 
solution exists to cure this problem: One just needs to flip a sign in one of the terms 
of the TeVeS model ITSl [T6l , and merely accept that photons can be slower than 
gravitons |29|. If all the field equations are ensured to remain hyperbolic, with a 
common time direction, no causal paradox can be caused by such a situation. 

The vector of action (fTsT l is assumed to be timelike in the TeVeS model, there- 
fore there a priori exists a preferred frame in which it takes the value C/^ = ( 1 , 0, 0, 0) 
II92II57I . References lfT5l[T6l argue that this can anyway be consistent with the high- 
precision tests of local Lorentz invariance of gravity in the solar system if this vec- 
tor is dynamical. However, Ref. |34| showed that the corresponding Hamiltonian 
is not bounded by below, precisely because of the kinetic term of this vector field. 
Therefore, the model is unstable. Other instabilities are also present in the slightly 
different model of Ref. |93J. 
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A final difficulty is related to the post-Newtonian tests of relativistic gravity, 
summarized in Sec. |4] above. In standard scalar-tensor theories Fig. Q shows 
that the linear matter-scalar coupling constant Oo should be small. In the TeVeS 
model, of the form ( fTSl l. the functions A(^) and B{(p) have been tuned to mimic the 
Schwarzschild metric of GR up to the first post-Newtonian order, even for a large 
matter-scalar coupling constant Oq. Therefore, the full plane (|ao|,j3o) of Fig. [T] 
seems now allowed by experimental data. However, binary-pulsar tests do not de- 
pend on A{(p) and B{(p) in the same way, and one actually gets basically the same 
constraints as in Figs. [TOl and [TTl notably the tight bounds on imposed by the 
pulsar-white dwarf binary PSR J1141 — 6545 of Fig.[TT] In conclusion, in spite of the 
tuning of A{(p) and B{(p) to mimic GR in the solar system, binary pulsars anyway 
impose that matter should be weakly coupled to the scalar field ll29l . This is quite 
problematic because we wish the same scalar field to give rise, at large distances, to 
the MOND acceleration ^/GMao/r whose magnitude is fixed. Figure [T2]illustrates 
schematically the difficulty (more precise discussions are given in Ref. ||29l ): If we 
wish the function / (((9^^)^) entering action ( fTSl ) to have a natural enough shape, 
either scalar-field effects are too large to be consistent with binary-pulsar tests, or 
the solar-system size is already large enough for all planets to be in the MOND 



Fig. 12 Derivative of the 
function / defining the kinetic 
term of the scalar field in ac- 
tion (ID, such that the MOND 
dynamics is predicted at large 
distances and Newtonian 
physics at small distances. In 
the first case, /' has a nat- 
ural shape and the MOND 
acceleration is obtained at the 
expected distance \/ GM/ao, 
but the scalar-field contri- 
bution to the acceleration 
is too large to be consistent 
with binary-pulsar tests. In 
the second case, /' has the 
same natural shape and the 
force caused by the scalar 
field is small enough when 
r — > 0, but all planets should 
undergo an extra MOND ac- 
celeration V CMao I r. In the 
third case, the scalar force 
is small enough in the solar 
system and the MOND law 
is predicted beyond the char- 
acteristic distance ^GM /ciq, 
but the shape of /' is ex- 
tremely unnatural. 
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regime and undergo a small ^GMgq / r acceleration in addition to the Newtonian 
GMjr^ . This would be ruled out by tests of Kepler's third law II103II . The only re- 
maining solution would be to impose a small enough a^GM/r^ acceleration caused 
by the scalar field within the solar system, and the MOND law ^GMoq / r at large 
distances, but this would correspond to the quite unnatural shape of the function 
/' displayed in the third panel of Fig. [12] Although this is not yet excluded ex- 
perimentally, it would suffice to improve binary-pulsar constraints by one order of 
magnitude to rule out this kind of fine-tuned model, because one would need a shape 
of function /' violating the consistency conditions (fT9l ). 



6.4 Nonmininal couplings 

Another possible way to modify gravity is inspired by the behavior of extended 
bodies within GR itself, as illustrated in Sec. |5] above. One may assume that matter 
is not only coupled to the metric but also nonminknally to its curvature. Let us thus 
consider an action of the form 



167rG* 



tter [matter; = fislv^Rl^vpXaRl^ivp^- • Oj > 

(20) 

which is Lorentz-invariant and satisfies the Einstein equivalence principle because 
all matter fields are coupled to the same tensor g^v (although it is built from the spin- 
2 metric g*^y and its curvature tensor). An immediate bonus of this class of theories 
is that they reduce to GR in vacuum, therefore standard solutions for the metric y 
remain valid, notably the Schwarzschild solution for spherically symmetric config- 
urations. The only difference with GR is that matter is no longer minimally coupled 
to g* y, and therefore that its motion can be changed. Reference ||29l showed that 
it is a priori possible to reproduce the MOND dynamics within this framework. 
However, the same theorem by Ostrogradski ISTlll 161 that we mentioned in Sec.|2] 
suffices to conclude that these models are generically unstable, and this can indeed 
be checked explicitly. 

This instability can fortunately be avoided by slightly complicating the model in 
vacuum. Instead of pure GR, let us assume a scalar-tensor theory in vacuum, with a 
negligible scalar mass nicp at galaxy scale. One may thus define an action 



c3 



S = J d^x^^ {R* -2g^''d^(pdy(p^2ml(p^} 

^'-'matter 

and show that specific functionals A[(p] = A{(p,dn(p) and B[(p] = B{(p,d^(p) allow 
us to reproduce both the MOND dynamics and the right amount of light deflection 
by galaxies. The generic Ostrogradskian instability is avoided because the matter 
action involves second time-derivatives of the scalar field only Hnearly ||29l . Such a 
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model is a priori as good as those which were previously proposed in the literature, 
and it is much simpler to analyze. Because of its simplicity, it has thus been possible 
to study its behavior within matter (an analysis which is usually too difficult to 
perform in other models). The bad news is that the scalar field equations do not 
remain hyperbolic within the dilute gas in outer regions of a galaxy. Therefore, this 
class of models is finally also ruled out to reproduce the MOND dynamics. 

On the other hand, this class (l2Tl i is able to reproduce the Pioneer anomaly (if 
confirmed) in a consistent relativistic theory of gravity, without spoiling the other 
well-tested predictions of GR. This anomaly is actually a simpler problem than 
galaxy rotation curves. Indeed, one of the difficulties of MOND models was to pre- 
dict a force ^/M. In the case of the two Pioneer spacecrafts, we do not know how 
their extra acceleration towards the Sun, 5a « 8.5 x lO^'^m-s^^, is related (or not) 
to the mass Mq of the Sun. There are thus actually several stable and well-posed 
models able to reproduce this effect, that were constructed in Ref. 1291 using the 
previous analyses of ll67l l68l l69l . One of the simplest models reads schematically 



where one needs to impose a? < 10^^ to pass solar-system and binary-pulsar tests, 
X sa {aGMQ/c^)^{5a/v^) « a^(10^'*m)^ to_^f the Pioneer anomaly, and where 
the scalar mass niff, needs to be negligibly small at solar-system scale. Actually, 
refinements are necessary to define correctly this model when (jO — > 0, but we refer 
to 1291 for this discussion. 



7 Conclusions 

In these lecture notes, we have underlined that the study of contrasting alternatives 
to GR is useful to understand better which features of the theory have been tested, 
and to suggest new possible tests. One may modify either inertia (i.e., the matter ac- 
tion S'matter), gravity (i.e., the action defining its dynamics, S'gravity). or consider non- 
minimal couplings to curvature. We stressed that the dynamics of gravity directly 
influences the motion of massive bodies, and this was illustrated by constructing the 
Fokker action. 

Our study of scalar-tensor theories of gravity exhibited a qualitative difference 
between solar-system experiments (weak fields) and binary-pulsar tests (strong 
fields). While the former tightly constrain the linear matter-scalar coupling con- 
stant, the latter are nonperturbatively sensitive to nonlinear couplings, and notably 
forbid that the quadratic matter-scalar coupling constant be large and negative. We 
mentioned that the no-hair theorem imposes that the motion of black holes in scalar- 
tensor theories is the same as in GR, contrary to other massive bodies. We also gave 




(22) 
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a simple dimensional argument showing that finite-size effects are much larger in 
alternative theories than in GR. 

Finally, we illustrated that the MOND phenomenology may a priori be repro- 
duced within any of the above classes of alternative theories (modified inertia, mod- 
ified gravity, couplings to curvature), but that all the proposed models present sev- 
eral experimental and theoretical difficulties. However, many possible routes remain 
possible, and generalizations of Einstein-aether theories |65 66 52l l64l[TT7l able 
to reproduce MOND deserve to be studied, notably from the point of view of sta- 
bility and well-posedness of the Cauchy problem. Of course, the simplest solution 
to account for galaxy rotation curves and cosmological data seems merely to accept 
the existence of dark matter. Instead of modifying gravity, one may then even de- 
vise some "modified dark matter" so that it reproduces the MOND successes while 
keeping the standard cold dark matter cosmological scenario II2TI l22l l23l . 
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